Electromagnetically induced transparency of an interacting cold Rydberg ensemble 
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We study electromagnetically induced transparency (EIT) of a weakly interacting cold Rydberg 
gas. We show that the onset of interactions is manifest as a depopulation of the Rydberg state and 
numerically model this effect by adding a density-dependent non-linear term to the optical Bloch 
equations. In the limit of a weak probe where the depopulation effect is negligible, we observe no 
evidence of interaction induced decoherence and obtain a narrow Rydberg dark resonance with a 
linewidth of <600 kHz, limited by the Rabi frequency of the coupling beam. 

PACS numbers: 03.67.Lx, 32.80.RrQ, 42.50.Gy 



Ensembles of Rydberg atoms display fascinating many- 
body behavior due to their strong interactions [l| . These 
interactions lead to interesting cooperative effects such 
as superradiance 0, 0, 0] and dipole blockade 0, H, 0, 11] , 
which may provide the basis for applications such as 
single-photon sources [1] and quantum gates [l^, [ll|, ■ 
For quantum information applications one is interested in 
the coherent evolution of the ensemble. Coherent excita- 
tion of Rydberg states has been achieved using adiabatic 
passage [13|, |lj]. Rabi oscillations between ground and 
Rydberg states with dipole-dipole interactions have been 
observed [l^, Also, the coherence of a Rydberg en- 
semble has been measured directly using a spin echo tech- 
nique (l7| . In most experiments on ultracold Rydberg 
gases, the Rydberg atoms are detected indirectly follow- 
ing field ionization and subsequent detection of electrons 
(or ions) using a micro channel plate (MCP). However, 
recently we demonstrated non-destructive optical detec- 
tion of Rydberg states in room temperature Rb vapor 
using EIT [23,IIl|. The same technique was sub- 
sequently used to detect Rydberg states in a Sr atomic 
beam |22i] . Rydberg EIT has number of potential appli- 
cations, for example, a Rydberg EIT medium displays a 
dc electro — optic effect many orders of magnitude larger 
than other systems 23] , and Rydberg EIT enables direct 
measurement of the coherence of the Rydberg ensemble. 

In this work we demonstrate EIT involving Rydberg 
states in an ultra-cold atomic sample. We show that in- 
teractions between Rydberg atoms lead to a rapid depop- 
ulation of the Rydberg state, and that the EIT spectra 
are extremely sensitive to small changes in the interaction 
strength. For example, changing the principal quantum 
number of the Rydberg state from n to n -t- 1 produces 
a significant change in the EIT spectrum. By using a 
weak probe, where the probability of populating the Ry- 
dberg state is low, we can eliminate this depopulation 
effect and obtain narrow Rydberg dark resonances with 
a linewidth of < 600 kHz, limited by the Rabi frequency 
of the coupling beam. 

The experimental setup and simplified level scheme are 
shown in figure 1 (a) and (b) respectively. A probe and 
coupling beam are combined using dichroic mirrors and 
counter-propagate through a cloud of laser-cooled ^^Rb 
atoms. The polarization of the beams are chosen to max- 




FIG. 1: (Color online) (a) Schematic of the experimental 
setup. A strong coupling beam (blue) stabilized to the 5p 
^P3/2 riD (or nS) transition counter-propagates a weak 
probe beam (red) through a cold sample of *®Rb released 
from a magneto-optical trap (MOT). The probe scans over 
a range of ±20 MHz about the _F = 3 — » F' = 4 resonance 
in 480 /is. (b) Simplified level scheme, Rabi frequencies and 
decay rates used to model the Rydberg EIT. Levels |4) rep- 
resent states other than state |2) to which the Rydberg state 
1 3) can decay. 



imize transition strengths {a'^-a'^). The probe beam is 
derived from a diode laser at 780.24 nm stabilized to 
the |1) ^ |2), 5s 2Si/2(f = 3) ^ 5p ^P^ (F' = 4) 
transition using polarization spectroscopy |24| . The fre- 
quency is controlled using two successive acousto-optic 
modulators (AOMs). The beam is then coupled into 
an optical fiber, the output of which is coUimated and 
passed through an aperture resulting in a beam of almost 
uniform intensity and radius of 0.75 mm. The coupling 
beam is derived from a commercial frequency doubled 
diode laser system (Toptica TA-SHG) at A=483-487 nm 
resonant with the |2) |3), 5p ^Pg/a {F' = 4) ^ nd 
^1)5/2 (or ns ^Si/2) {F") transition. The frequency of 
the coupling laser is stabilized using a Rydberg EIT sig- 
nal produced in a magnetically shielded rubidium vapor 
cell at room temperature [l^l ■ 
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We load approximately 10^ atoms in a magneto-optical 
trap (MOT) in 5 s. A subsequent 10 ms optical molasses 
phase cools the atoms to a temperature of approximately 
20 /iK and density 10^" cm^^, following which the cooling 
light is extinguished. An optical pumping pulse of vari- 
able duration prepares the atoms in the F = 3, mp = 3 
state. After 3 ms free expansion the pump and probe 
beams are turned on and the frequency of the probe is 
swept across the D2 F ^ 3 ^ F' — 4 transition. The 
sweep covers ± 20 MHz about the resonance in 480 /iS. 
Absorption measurements show that there is negligible 
density decrease over the duration of the scan and neg- 
ligible optical pumping into the lower [F = 2) hyperfine 
ground state for the probe powers of interest. The probe 
beam is incident on a fast photodiode which records the 
transmission through the sample. This technique of mea- 
suring the EIT spectra from a cold atomic sample within 
a single experimental cycle offers a fast repetition rate. 

In figure 2 we show a series of EIT spectra for a fixed 
probe power of 3.6 /zW corresponding to 0.13 /sat, where 
/sat is the saturation intensity on the |1) — > |2) transition. 
The absorption data has been corrected for the measured 
change in probe intensity through the optical fiber as a 
function of AOM frequency. The effect of varying atom 
density, p, and principal quantum number, n, are shown 
in figures 2 (a) and (b), respectively. The atom den- 
sity, p, is changed by heating the atom cloud (of initial 
density po) with the probe beam prior to the EIT scan. 
In figure 2 (b) the power of the coupling laser is varied 
to make the coupling Rabi frequency the same for each 
scan. For low p or low n a characteristic EIT spectrum 
is observed with a slight asymmetry. As p or n are in- 
creased this asymmetry becomes more pronounced until 
the transparency peak is lost completely. The lineshape 
for n > 22 does not change dramatically with n as il- 
lustrated by the similarity to the n — 26 and n = 22 
data sets, the front curves in figure 2 (a) and (b), respec- 
tively. The data presented in figure 2 is suggestive of a 
density-dependent depopulation of the Rydberg state. 

By varying the probe power we can change the frac- 
tional probability that an atom is excited to the Rydberg 
state. In figure 3 we show the transmission spectrum 
for the |3) = 19d state for two probe powers of 450 nW 
(0.016 /sat) and 3.6 /xW (0.13 /sat)- In this case we 
present double scans where the probe laser detuning is 
first increased and then decreased through the resonance 
in a total time of 960 ps. We observe that at the higher 
probe power there is significantly less absorption on the 
second scan than the first indicating that atoms excited 
to the Rydberg state are being lost from the EIT sys- 
tem. For these parameters only a small fraction of the 
population, a few percent, is transferred into and out of 
the Rydberg state for a time of order 10 ps. However, to 
model the data (see next paragraph) we find that >70% 
of the atoms are lost from the 3-level EIT system. This 
substantial loss cannot be explained by spontaneous de- 
cay from Rydberg state (since the spontaneous lifetime 
is also of order ~10 pa ^a] and decreases with larger n) 
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FIG. 2: (Color online) (a) Transmission of a 3.6 pW probe 
beam as a function of probe laser detuning for the *^Rb 
|1) 1 2) —> 1 3) — 26d system for varying atomic density, 
(b) EIT spectra for n=19-22. The Rabi frequency of the 
|2) j3) transition is kept the same for each state. In both 
plots the onset of density-dependent behavior is evident in 
the emergence of an EIT hneshape at lower values of n and p 



or by optical pumping to the lower hyperfine state. We 
therefore deduce that there is an interaction induced loss 
from the system happening on a much faster time scale 
than the Rydberg state lifetime. 

To describe this interaction-induced depopulation we 
add non-linear density dependent decay and decoherence 
terms to the standard three-level optical Bloch equations 




iha — [//, a] — 7(7 (1) 

where a is the density matrix, H the atom-light interac- 
tion Hamiltonian and 717 the decoherence matrix. The 
depopulation of the Rybderg state is modeled by intro- 
ducing a fourth level |4) to represent the reservoir of 
states other than |2) to which an atom in Rydberg state 
1 3) could decay, as shown in figure 1 (b). The decay of 
the Rydberg state |3) is given by 730-33-1-73(713 where the 
first term is due to spontaneous emission and the second 
term is a density-dependent interaction induced loss to 
the reservoir states |4). The spontaneous decay rates of 
states 1 2) and |3) are given by 72 and 73, respectively. We 
assume that 50% of the spontaneous decay from |3) pop- 
ulates state 1 2) and 50% is lost to the reservoir |4), how- 
ever this weighting makes little difference to the lineshape 
as it is the non-linear term, 73, that dominates. Atoms 
decay from |4) to |1) at a rate 74. The coherences are as- 
sumed to decay at half the rate of the populations. The 
decoherence rate between |1) and |3), 731, is included as 
a fitting parameter which depends on the combined laser 
linewidth and dephasing due to Rydberg atom interac- 
tions. The optical Bloch equations are then solved using 
a time varying probe detuning to match the experimental 
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sequence. As the probe is swept through the two-photon 
resonance the Rydberg state is first populated and then 
depopulated by chirped adiabatic passage, in contrast to 
conventional adiabatic passage schemes where the trans- 
fer is unidirectional p^ll4| . The predictions of the model 
are compared to the experimental spectra in figure 3. 
The model accurately reproduces the observed spectra 
for different values of the probe laser intensity (all other 
parameters remain the same). When the probe intensity 
is increased the population of the Rydberg state (sohd 
lines in the lower panel of figure 3) also increases and 
there is a faster decay into the reservoir states (dotted 
line in the lower panel of figure 3). For a higher probe 
power of 3.6 /xW the model predicts that about 3% of the 
population is transferred into the Rydberg state and after 
the double scan sequence almost 80% of the population 
is in other states, |4). 
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FIG. 3: (Color online) Top: Double scan BIT for the jl) 

1 2) — * 1 3} = 19d system for probe powers of 450 nW (red) 
and 3.6 fiW (green). The model fits are shown as the solid 
black lines. Bottom: The populations of the Rydberg states 

1 3) x25 (solid) and reservoir states |4) (dashed) as a function 
of time during the scan 

There are a number of surprising features of this 
density-dependent depopulation effect. First, we ob- 
serve similar decay rates for nearby s-states. Second, 
the observed decay is relatively fast for such low Rydberg 
atom densities of order 10® cm~^. For example, the time 
scale associated with near-resonant dipole-dipole energy 
transfer processes [26| at this density would be 10 to 20 
times slower. Also, our simple non-linear optical Bloch 
equation model begins to fail for larger n. For example, 
for the high density 26d trace shown in figure 2 (a) no 
evidence of FIT is observable whereas the model predicts 
a pronounced shoulder after the laser has passed through 
the two-photon resonance. One possible explanation is 
that the decay out of the KIT system is driven by su- 
perradiant decay of the Rydberg state to nearby f and 
p states. The potential importance of superradiance in 
frozen Rydberg gases particularly involving lower n states 
has been highlighted in recent work 0, [23| ■ In general su- 
perradiance is more significant for low n as the decay rate 



to nearby Rydberg levels scales as 1/n^, however there 
is also a geometrical effect which depends on the wave- 
length of the superradiant transition. As, in comparison 
to the size of the cloud, D [l^l- For the range of n studied 
in this work {n — 19 to 26) As is less than D for decay 
to the closest f states i.e. the atoms are spatially spread 
over a range greater than the transition wavelength such 
that cooperative effects are small. However, As is of order 
or greater than D for decay to the closest p states i.e. all 
atoms are within one wavelength and the emission rate 
is greatly enhanced. In this critical region the geometric 
cooperativity parameter [27t increases rapidly with in- 
creasing n which could explain the dramatic evolution of 
EIT spectra shown in figure 2 (b) . The dependence of su- 
perradiance on As and D is supported by the similarity of 
figure 2 (a) and (b) where for both data sets the value of 
XsD varies by a factor of 1.6 between the front and back 
trace psj . Finally, we note that one can obtain a better 
fit to the high density and higher n data by coupling the 
optical Bloch equations to an equation describing the su- 
perradiance field energy density. This model reproduces 
the lineshapes observed in the experiment, in particular 
the complete loss of the EIT shoulder, however as the 
time scale of the decay of the superradiance field is not 
known, this model does not provide useful quantitative 
information. 

Lastly, we focus on the weak probe limit where the pop- 
ulation of the Rydberg state becomes negligible. This 
limit highlights the power of EIT as a non-destructive 
probe of Rydberg gases. Figure 4 shows the transmission 
spectrum for the |3)=26d state using a low probe power 
of 200 nW (0.007 /sat)- In this case no loss or additional 
decoherence of the Rydberg state is observed between 
the forward and back scan. We obtain a typical EIT 
linewidth of 600 kHz. The EIT linewidth is made up of 
two terms; a term proportional to the coupling Rabi fre- 
quency squared, f2|2 (i-e- the Autler-Townes splitting), 
and a term corresponding to the coherence between the 
ground state and the Rydberg state, 731. This coherence 
term is the sum of the combined laser linewidth and inter- 
action dependent dephasing mechanisms. From the fit we 
obtain n32/2n = 1.8 ± 0.3 MHz and 731 = 100 ± 20 kHz, 
i.e. for our parameter the linewidth is dominated by the 
Autler-Townes term. We observe no change in 731 over 
the range n = 19 to n = 26 showing that the combined 
laser linewidth dominates over any interaction induced 
dephasing of the Rydberg state. Note that such long co- 
herence times cannot be observed in non-Rydberg cas- 
cade EIT due to the short lifetime of the upper state. 
The frequency resolution obtained using this method is 
comparable to the best obtained on single atoms [lH| and 
considerably better than the resolution obtained for other 
ensembles experiments [16, . 29., 30] , where interaction ef- 
fects tend to broaden the linewidth to greater than one 
MHz. This narrow resonance demonstrates a significant 
advantage of EIT for precision spectroscopy or electrom- 
etry, i.e. that the observed transparency is sensitive to 
the Rydberg energy level without transferring significant 
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population into it. From the model we find that the 
probability for an atom to be in the Rydberg state to be 
~0.6% giving p < 10'' cm~'^ which corresponds to a total 
Rydberg population of less than 10^ atoms. 
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FIG. 4: (Color online) Top: Double scan EIT for the |1) 
1 2) —> 1 3) = 26d system. The solid line (black) is the the- 
oretical line of best fit giving an EIT linewidth of 0.58 ± 
0.04 MHz. Bottom: Population of the Rydberg state as a 
function of time for the duration of the scan. 



In conclusion, we have shown that EIT is a useful tech- 
nique to study the properties of interacting Rydberg en- 
sembles. We find that interaction effects are manifest as 
a rapid depopulation of the Rydberg state. This den- 
sity dependent loss can be explained using a simple non- 
linear optical Bloch equation model. The character of 
the loss is similar for changes in either the atomic den- 
sity or the principal quantum n of the Rydberg state, 
suggesting a superradiant decay mechanism. In the limit 
of weak probe power, where the depopulation effect is 
negligible we observe no evidence of interaction induced 
decoherence and the linewidth is limited by the Autler- 
Townes splitting. In this limit EIT provides a sensitive 
non-destructive probe of a cold Rydberg ensemble allow- 
ing high-resolution spectroscopy of Rydberg states. 
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